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ABSTRACT 

According to TV-body numerical simulations, the radial pairwise velocities of galaxies 
C^Q ' have negative skewness in the quasi-nonlinear regime. To understand its origin, we 

calculate the probability distribution function of the radial pairwise velocity using the 
Zel'dovich approximation, i.e., an analytical approximation for gravitational cluster- 
ing. The calculated probability distribution function is in good agreement with the 
result of TV-body simulations. Thus the negative skewness originates in relative mo- 
tions of galaxies in the clustering process that the infall dominates over the expansion. 

Key words: cosmology: theory — large-scale structure of universe 

P ' 1 INTRODUCTION 

CZ3 . Let us consider two galaxies separated from each other by the vector r. If their peculiar velocities are u\ and 1x2, the pairwise 
' velocity v is defined as 

^ ! 

= 1x2 — «i = v«n« + v±n±. (1) 



Here the unit vectors toii and n± are parallel and perpendicular to the separation vector r. We study probability distribution 
' functions (PDFs) of vu and v±, which are respectively denned as radial and transverse components of the pairwise velocity 
(— 00 < v« < +00, ^ v± < +00). 

The radial pairwise velocity is of fundamental importance in observational cosmology (Peebles 1993, section 20). It serves 
as a probe of the matter distribution. Its PDF is essential to converting observed data in the redshift space into those in the 
real space. 

For the quasi-nonlinear regime of the present-day universe (r = |r| ~ 10 1 -10 2 Mpc), TV-body numerical simulations have 
shown that the radial-velocity PDF exhibits a negative average and negatively skewed pronounced tails (Efstathiou et al. 
1988; Fisher et al. 1994; Zurek et al. 1994; Magira, Jing & Suto 2000). These properties are attributable to coherent motions 
of galaxies approaching each other, since gravitational clustering is the predominant process in the quasi-nonlinear regime. 

However, with the results of TV-body simulations alone, we cannot understand how the negative average, negative skewness 
and pronounced tails of the radial- velocity PDF are related with gravitational clustering. To clarify this relation, an analytical 
treatment is desirable. 

The linear perturbation, i.e., the most convenient analytic tool, in insufficient to study the above-mentioned properties 
of the radial-velocity PDF. This is because these properties are nonlinear. The linear perturbation merely amplifies peculiar 
velocities of all the galaxies by the same factor. Thus the radial-velocity PDF retains its Gaussianity in the initial stage, where 
motions of the individual galaxies are random and independent. Although the leading terms of the average (v«) and standard 
deviation («|) happen to be squares of the linear perturbation, those of the higher-order moments, e.g., (hm) and are 
nonlinear perturbations (see Juszkiewicz, Fisher & Szapudi 1998). 

Juszkiewicz et al. (1998) used a second-order Eulerian perturbation and reproduced successfully the moments up to the 
third order, (v\\), (v?\), and as well as the negatively skewed PDF. However, the PDF was obtained with a conventional 
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but unjustified ansatz, instead of the rigorous theory that requires complicated calculations. Seto & Yokoyama (1998) used a 
first-order Lagrangian perturbation, i.e., the Zel'dovich approximation (ZA; Zel'dovich 1970). Although their radial- velocity 
PDF exhibits pronounced tails, they failed to reproduce the negative average and had to shift their PDF by hand in order 
to achieve agreement with A^-body simulations. Overall, the mechanism that determines the radial-velocity PDF is still 
controversial. 

We use ZA to reexamine the PDF. It would be surprising if ZA really reproduced the pronounced tails but not the negative 
average The former is associated with the higher moments such as (v*)- For the quasi-nonlinear regime, ZA is a good 

approximation (Yoshisato, Matsubara & Morikawa 1998). Each of the galaxies is assumed to move along the gravitational 
force that is determined by the initial density distribution. Thus the galaxies cluster together. The velocity and density fields 
undergo a nonlinear evolution, which is expected to cause the negative average, negative skewness and pronounced tails of 
the radial-velocity PDF. 

We have to be careful about the extent to which ZA can describe gravitational clustering. Of particular importance is 
the range of the separation r where ZA is applicable. It has been known that ZA is applicable only to a limited range, but 
the applicable range has not been estimated so far. The quantitative estimation is done for the first time in this paper. We 
subsequently demonstrate that Seto & Yokoyama (1998) made their ZA analysis outside its applicable range. 

Throughout the main text, we adopt the Einstein-de Sitter universe with a Hubble constant Ho = 50 km s _1 Mpc^ 1 for 
comparison with the previous works. In Section 2, an analytic form of the joint PDF for the radial and transverse velocities is 
derived with ZA. The applicable separation range of ZA is estimated in Section 3. We compare ZA with A r -body simulations 
in Section 4. The origin of the negative average and negative skewness is discussed in Section 5. We discuss the existing 
relevant models in Section 6. The conclusion is presented in Section 7. The application to other cosmological models and so 
on are discussed in Appendices. 



2 PAIRWISE VELOCITY DISTRIBUTION 
2.1 Analytical form 

Using ZA, we newly derive an analytic form of the joint PDF for the radial and transverse velocities at the separation r 
and the time t, P(vu,v±\r,t). It is assumed that the vectors v\\n\\ and Diiii have isotropic distributions. The joint PDF is 
normalized as 

P (v\\,v±\r,t) dv\\2Trv±dv± = £(r,t) + 1, (2) 

where £(r,t) is the two-point correlation function for the number density of galaxies. The radial- velocity PDF, P(v\\\r,t), is 
obtained by integrating the joint PDF over the transverse velocity v± . 

Let us consider a relative motion of two galaxies. Their separation vector r — r(t) and pairwise velocity v = v(t) are 
written with ZA as 

D 

r = x-2 - xi = n + —Vi, 

Di b (3) 

V — U'2 — Ml = Vll nil + V±Tl± = —r-Vi. 

" " Di 

Here D = D(t) oc t 2 ^ 3 is the linear growth factor of the density fluctuation. We set D = 1 for the present-day universe. The 
suffix i indicates quantities at the initial time t = ti. Since the galaxies move along straight lines, the relative motion lies on 
the plane determined by the vectors n and Vi. 

With an increase of the time from t to t + dt, a galaxy pair is defined to move from (v\\,Vj_,r) to Mi,v'j_,r'). From 
conservation of the number of galaxies, we have 

P (vy , v'±\r', t + dt) dv\\2irv'j_dv'j_<iirr' 2 dr = P (v\, v±\r, t) dv\\2nv±dv±4:iTr 2 dr. (4) 
From equation the time evolution of dvu 2nv±dv±4nr 2 dr is obtained as 

dv'\\ 2nv' ± dv' ± 4TTr' 2 dr' = ^1 + S^dt^J dv^2nv±dvj_4nr 2 dr. (5) 
The solution is 



dv^2-KVj_dv±4:TTr dr = ( ) dv^ >i 2Trvj_ ti dv± t i4:iTr i dri. (6) 
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Thus the joint PDF at the time t is related to the PDF at the initial time ti as 

P(v h v±\r,t) = P(t)||,i,t;±,i|r 4 ,ti). (7) 
From equation the initial separation n is obtained as 

= \/f r - +f5^1 • ( 8 ) 



Likewise, the initial pairwise velocities Wiij and v± ti are obtained as 
Di frv\\ Dvf { +v 2 ± \ 

)' (9) 



Vti = 

D n 



D \ ri D 



(see also Appendix A). Seto & Yokoyama (1998) used ZA to derive directly an analytic form of the radial-velocity PDF, 
P(u|l \r, i). We prefer the joint PDF, P(v\\ , v± \r, t), which has a more simplified form and provides us with more straightforward 
information about, e.g., the origin of the negative skewness in the radial velocity (Section 5). 



2.2 Initial condition 



Here we determine the initial condition for P(v\\,v±\r,t). Although the same formulation was adopted in Seto & Yokoyama 
(1998), we give a more complete description for our discussion in Sections 3 and 5. The basic assumption is that the initial 
density fluctuation is random Gaussian while the initial peculiar-velocity field is homogeneous, isotropic, random and vorticity- 
frce. 

Since the initial density fluctuation is random Gaussian, the corresponding PDF of the radial and transverse velocities is 
also Gaussian (Fisher 1995), 



P (*y||,i,«x,i|n,tj) 



(27r)3a2(r 1 )al(r l ) 



: exp 



+ 



2 V ffjf(ri) ffi(ri) 



(10) 



Thus we only have to obtain the dispersions of the initial pairwise velocities o-y (n) and a±(ri). The averages and (w±,i) 

do not have to be considered because they vanish in the limit Di — > 0. 

Since the initial peculiar-velocity field is homogeneous, the pairwise-velocity dispersions are derived from the one-point 
velocity dispersion and the two-point velocity correlation: 



(11) 



0-jj(r») = 2{u| ]i (xj)) a - 2(tt|| li (xi)ii|| ii (x< + rj)) s 

a ±( r i) = ^{ U \,i{ X i))s - 2(u ± ,i(Xi)u ± ,i(Xi + n)) a . 

Here (-) s denotes a spatial average and should not be confused with (■) that denotes an average over galaxy pairs; wm^ and 
u±,i are the radial and transverse components of the peculiar velocity Ui, respectively. We relate the velocity correlations to 
the power spectrum of the density fluctuation as follows (Gorski 1988). Using a theory for a vector field that is homogeneous, 
isotropic, random and vorticity-free (Monin & Yaglom 1975, section 12), we have 



(u|| li (xi)«|| )j (a!j + r,)) s 



jo (fen 

kiTi 



kiTi 
&i (J^i ) dki . 



Here jo and ji are first-kind spherical Bessel functions, 

sin x cos x 



. , , sin a; 

jo(x) = and ji{x) 



and £i(ki) is the power spectrum of the initial peculiar-velocity field, 



£i(ki)dki = - (\ui(x t )\ 2 ^ 



(12) 



(13) 



(14) 



Using the linear perturbation, we relate the power spectrum of the initial peculiar-velocity field £i(ki) to that of the initial 
density fluctuation Pj(fcj). The linear peculiar- velocity field is described by the linear density field as 



.Di 



u l (x i )=i— I -^5i(ki)exp(ikiXi) 



dki 



(2tt) 3 / 2 



(15) 
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where Si(ki) is the Fourier transform of the density contrast Si(xi) (Peebles 1993, section 21). The power spectrum Vi(ki) of 
the initial density fluctuation is 



(8i[ki)W)*) = v *(fa) ( 2?r ) 3 6 ( k * - fe 



(16) 



where * denotes the complex conjugate. From equations 114t - 116ll . the power spectrum £i(ki) of the initial peculiar- velocity 
field is obtained as 



SAl.;) [ ^ 



Vi{ki). 



(17) 



Therefore, the radial- and transverse-velocity dispersions at the initial stage are 



a \\( r ') = — ( TT 



8tt 
3 

3 



■ 3 jo (far i) + 6 
ji{fari 



ji {fari 



Vi(ki)dki = 



v-.a-.xik, = ( ^ 



RKn), 



(18) 



Rl(n). 



The initial condition for P(vu,Vx\r, i) is determined by equations HUH and 11811 . 

Since the peculiar velocity depends on the gravitational potential produced by the density fluctuation, equation 1181 
includes the information that galaxies are about to cluster together. In fact, the radial velocity dispersion a 2 (n) reflects the 
two-point velocity correlations {u\\,i(xi)u\\ t i(xi +ri))„. At a small separation n where galaxies are about to move in the same 
direction, the correlation is positive. At a moderately large separation where galaxies are about to cluster from the opposite 
sides, the correlation is negative. At a very large separation, the correlation is absent. 



2.3 Power spectrum of initial density fluctuation 

The power spectrum of the initial density fluctuation Vi(fa) is adopted from the standard cold dark matter model (Peebles 
1993, section 25): 

v + h) = r — f ^Sr (19) 

where a = 25.6 Mpc, /3 = 12 Mpc, 7 = 6.8 Mpc and v = 1.13 for the Einstein-de Sitter universe with Ho = 50 km s _1 Mpc - 
(Bond & Efstathiou 1984; Efstathiou, Bond & White 1992). The normalization factor B is determined from the temperature 
fluctuation of the cosmic microwave background: 

B - 5^ 4 l^rJ ■ (20) 

Here the quadrupole fluctuation amplitude Q r ms is set to be 9.5 fiK with the present-day temperature To = 2.73 K. The same 
initial power spectrum was adopted in JV-body simulations of Fisher et al. (1994) and Zurek et al. (1994), which are to be 
compared with our calculation, and also in the ZA calculation of Seto & Yokoyama (1998). Other initial power spectra and 
cosmological models are studied in Appendix B. 

Note that there are relations vnj oc Di, Vj_,i cc Di, o"i|(rj) oc Di, a±(ri) oc Di and Vi(fa) oc D\ in equations @, I18H and 
l|19jl . Thus the initial condition for P (v\\ , v± \r, does not depend on the values of Di and Di. 

The initial spectrum Viiki) is used without any modification. Coles, Melott & Shandarin (1993) suggested that ZA's 
performances are improved if small-scale fluctuations are removed from the initial spectrum by using a window function. 
Although this appears to be true in the calculation of the matter distribution, the amplitudes of pairwise velocity dispersions 
would be incorrectly underestimated as shown in Seto & Yokoyama (1998). Therefore, we do not adopt this modification. 



3 APPLICABLE RANGE OF ZEL'DOVICH APPROXIMATION 

Galaxies in ZA continue to move straight, and pass away each other even if they once cluster together. Since actual clustering 
galaxies are bounded gravitationally, ZA is no longer valid after this so-called shell crossing. At a fixed time, ZA is not 
applicable to small separations where the shell crossing is statistically significant. In Section 3.1, we analytically estimate the 
separation range where ZA is applicable. In Section 3.2, we make a numerical estimation on the basis of time evolution of the 
two-point correlation function. The applicability in the velocity range is also discussed in Section 3.3. 
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Figure 1. Graphical representation of equation 1221 at D = 0.1, 0.2, 0.4, 1, 2 and 10 as a function of the separation r. The dashed line 
denotes the left-hand side of the equation. The solid lines denote the right-hand side. They are in units of Mpc. ZA is applicable over 
the separations where the dashed line is well above the solid line. 



3.1 Analytical estimation 

Let us consider a pair of galaxies that are approaching each other. If its initial separation is r;, its typical initial radial velocity 
is estimated as Viu — —<r«(ri) — ~DiD~ 1 R^(n). From equation we have 

ri-r~--J.wii.i~ £-R\\(n). (21) 



The shell crossing occurs if r — 0. Thus ZA is valid as far as n 3> DD^ 1 Ru (n). Since this condition implies r ~ n, we replace 
n with r and obtain 
D 



r > jrR\\(r 



(22) 



This is ZA's applicable range. Fig. 1 compares the right-hand side of equation 1221 with the left-hand side as a function of 
the separation r at D — 0.1, 0.2, 0.4, 1, 2 and 10. ZA is applicable over the separations where the dashed line (the left-hand 
side) is well above the solid line (the right-hand side). It is evident that the applicable range is progressively narrow with a 
progress of the time. At D — 1, ZA is applicable to r 3> 10 Mpc. 



3.2 Numerical estimation 

The shell crossing affects the time evolution of the two-point correlation function. This is because, at a given separation, 
galaxies become less clustered after the onset of the shell crossing. Using ZA, we numerically calculate the two-point correlation 
function £ and its time derivative d^/dt at D = 1 (see equation The results are plotted as a function of the separation r 
in Fig 2. The time growth d^/dt of the two-point correlation function is maximal at r ~ 10 Mpc. This separation coincides 
with that for r = DD~ 1 R^(r) (fig. 1 in Section 3.1). Thus ZA's applicable range is r > 10 Mpc at D = 1. Any discussion on 
ZA has to be restricted within this range. 



3.3 Applicability in velocity range 

Even if the separation r is within ZA's applicable range, ZA is valid only for the velocity range 

v\\ « §r. (23) 

For U|| > _DD _1 r > 0, the initial radial velocity is negative as shown in equation ©. Such a galaxy pair has once clustered 
together with a negative radial velocity, experienced the shell crossing, and then turned to have a positive radial velocity. In 
the limit — » oo, equations (171- 1101 1 and 11811 yield a Gaussian PDF, 
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Figure 2. Two-point correlation function £ (a) and its time derivative d^/dt (b) at D = 1 in ZA. We show d^/dt in units of t 1 , where 
to is the present age of the universe. The abscissa is the separation r in units of Mpc. 



P (y\\,v±\r,t) 



exp 



2 i 2 



2R C 



where 



Stt 



This is overestimation for the positive radial velocity because of the shell crossing as stated above. 



(24) 



(25) 



4 COMPARISON WITH TV-BODY SIMULATIONS 

The radial-velocity PDFs are shown in Fig. 3 for r = 1, 5, 10 and 23 Mpc at D = 0.1, 0.2, 0.4, 1, 2 and 10 (solid lines). They 
are compared with the results of TV-body simulations of Zurek et al. (1994) and Fisher et al. (1994) at D = 1 (dashed lines). 
For reference, we also show Gaussian PDFs (dotted lines). 

Fisher et al. (1994) presented the radial-velocity PDF at D = 1 for £ = 0.1, which corresponds to r ~ 23 Mpc (fig. 2a). 
This separation is within ZA's applicable range. The PDF for £ = 0.1 in the A?-body simulation is in agreement with the PDF 
at r = 23 Mpc in ZA (fig. 3) . Both the PDFs exhibit a negative average as well as negative skewness. They are surely due to 
gravitational clustering of galaxies. 

The average (v\\), standard deviation ((v\\ — (w||)) 2 ) 1//2 and skewness ((wm — (vu )) 3 )/((«|| — ( v \\)) 2 ) 3 ^ 2 of the radial velocities 
are —190 km s _1 , 430 km s" 1 and —0.69, respectively, in the A^-body simulation. They are -133 km s"\ 462 km s" 1 and 
—0.36, respectively, in ZA. The velocity limit DD~ 1 r is 1150 km s . The average and standard deviation in the A'-body 
simulation are nearly the same as those in ZA. However, the skewness in the A'-body simulation is somewhat different from 
that in ZA. The radial- velocity PDF in the A'-body simulation has a more pronounced tail in the negative side. This is because 
gravitational acceleration of clustering galaxies is more significant than that assumed in ZA, which simply extrapolates the 
initial velocity Vi (equation |3J. 

For the separation r = 23 Mpc, the shell crossing becomes significant and ZA becomes invalid at D ~ 2 (fig. 1). The PDF 
becomes symmetric (fig. 3). Eventually, almost all the galaxy pairs experience the shell crossing. The PDF becomes Gaussian 
as observed at D = 10. 

Zurek et al. (1994) presented the radial-velocity PDF at D — 1 for r ~ 1, 5 and 10 Mpc. These separations are outside 
ZA's applicable range at D = 1 (figs 1 and 2). Fig. 3 shows the evolution of the PDF. At D = 0.1, ZA is valid. The PDF 
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Figure 3. Radial- velocity PDFs. Solid lines: ZA calculation at D = 0.1, 0.2, 0.4, 1, 2 and 10 for r = 1, 5, 10 and 23 Mpc. Dashed lines: 
JV-body simulations at D = 1 of Zurek et al. (1994) for r = 1-2, 5-6 and 10-11 Mpc and of Fisher et al. (1994) for r = 23 Mpc. Dotted 
lines: Gaussian PDFs. The abscissa is in units of the radial-velocity dispersion <T||(r) in ZA, which is 462 km s _1 at D = 1 for r = 23 
Mpc. 



exhibits a negative average and negative skewness. With increasing D, the shell crossing becomes significant, ZA becomes 
invalid, and the PDF becomes symmetric. At D = 1, the PDF in ZA is quite different from that in the A-body simulation. 
While ZA yields (v\\) = +8, +9 and — 54 km s for r = 1, 5 and 10 Mpc, respectively, the A-body simulation yields 
= -280, -355 and -276 km s _1 . 

Seto & Yokoyama (1998) made a ZA calculation for r = 1-11 Mpc at D = 1. Since ZA and A-body simulations yield 
different radial-velocity PDFs, they concluded that gravitational clustering is not incorporated completely in ZA. Since the 
PDFs for r ~ 1 and 5 Mpc in ZA and A-body simulations exhibit exponential tails, In P(v\\ \r, t) oc — |t>|| |, they also concluded 
that simple kinematics as ZA results in these tails. However, the separations r = 1-11 Mpc at D — 1 are outside ZA's applicable 
range. The exponential tails at r = 1 and 5 Mpc in ZA are spurious and attributable to galaxies that have experienced the 
shell crossing for the positive side of the PDF and to galaxies that are still in gravitational infall for the negative side (see 
also Section 6). 



5 ORIGIN OF NEGATIVE SKEWNESS 

The negative average and negative skewness of the radial-velocity PDF originates in galaxy clustering, i.e., the infall vv, < 
dominates over the expansion vn > 0. We discuss how this process is incorporated in ZA. Since ZA is valid at DD~ X (r) <C r, 
we separately study the velocity ranges |wm| <C DD~ Rn, DD~^R\\ <C |t>ii| <m DD~ 1 r and DD~ 1 r <ti \vn\, where DD~ _Ry = 
495 km s _1 and i>D~ x r = 1150 km s' 1 at D = 1 for r = 23 Mpc. 

First, we study the velocity range DD^R^ < |«||| < DD~ 1 r. Since most of the galaxy pairs have relatively small 
transverse velocities vj_ ~ DDT 1 a^{r) ~ _DD~ 1 i?x(r) ~ DD^ 1 R\\(r), equations © and © are simplified as 

n * r(l-^), (26) 
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Figure 4. (a) Dependence of the initial pairwise- velocity dispersions an (ri) (solid line) and a±_ (r.;) (dashed line) on the initial separation 
ri. The dispersions are in arbitrary units. The separation is in units of Mpc. (b) Two Gaussian PDFs with different standard deviations. 



Dt , Di 

Vn t ~ —v\\ and v±,i ~ — v±. 
ii, ^ ii D 

Thus the initial separation depends on the sign of the present radial velocity v\\ 

n(v\\ = +V) < n(v\\ = -V), 



(27) 



(28) 



where V > 0. For the initial separation r*, the radial-velocity dispersion ctii (t,) has the following dependence as shown in Fig. 
4a: 



ay (r-) < (Tii (r") for r • < r" < 100 Mpc. 



(29) 



The initial PDF, f^viij, v± t i\ri,ti), is Gaussian with the standard deviation a^(ri). Equation 1271 implies \vu i\ 2> crii^i). For 
such an argument that exceeds the standard deviation, a Gaussian function with a larger standard deviation yields a larger 
value as shown in Fig. 4b. Hence we expect 



P (v\\,i,vj.,i\r'i,u) < P (v\\,i,vj_ t i\r",ti) for r- < r" , 
which leads to 

P (wy = + V, v±\r,t) < P («n = - V, vx \r, t) . 



(30) 



(31) 



Thus the infall Vu < dominates over the expansion vu > in this velocity range. The PDF is accordingly asymmetric. 

Second, we study the velocity range |«ii| <C DD~ 1 R^. Since the absolute value of vu tends to be much smaller than v±, 
the sign of v« is not important to the values of ri and v\\j in equations ||S) and ©. We expect P(vu = +V,v±\r,t) ~ P(v\\ = 
— V, v±\r, t). The radial velocities among galaxies with small peculiar velocities retain the initial Gaussian character. These 
galaxies have not moved significantly from their initial positions. 

Third, we study the velocity range |«ii| S> DD~ 1 r. Equations © and @ yield ri ~ DD^ 1 ^^] and hence v\\ i — 
— 73i_D _1 |w|| |, which in turn yields P(v\\ = +V,v±\r,t) ~ P(v\\ = —V,V±\r,t). This is due to shell crossing as discussed in 
Section 3.3. For the real universe, the radial-velocity PDF is expected to be asymmetric in this velocity range. 

Therefore, the negative average and negative skewness in ZA stem from the inequality P(vn = +V, v± \r, t) < P(i)\\ = 
— V, vj_ \r, t) in the velocity range DD' 1 ^ <^ \vn \ DD~ 1 r. The existence of this velocity range is assured by the condition 
for ZA's applicability, DD' 1 ^ <g r. Whenever ZA is applicable, its radial-velocity PDF has a negative average and negative 
skewness. 
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The dependence of the radial-velocity dispersion cr«(ri) on the separation ri in equation 12911 is also an assured behavior. 
As discussed in Section 2.2, the radial- velocity dispersion cry (r^) reflects the difference between the one-point velocity dispersion 
(tt| i{xi)) s and the two-point velocity correlation (u\\ t i(Xi)u\\ t i(xi + i"j)) s . They are identical at n = 0. With an increase of 
the separation n, the two-point correlation decreases slowly in accordance with motions of galaxies that are about to move 
coherently and cluster together. 

6 DISCUSSION ON RELEVANT MODELS 

The pairwise velocity was studied by Juszkiewicz et al. (1998) using a second-order Eulerian perturbation. At the separation 
of 20.8 Mpc, the average, standard deviation and skewness of the radial velocities are -200 km s"\ 430 km s" 1 and -1.00, 
respectively. They are —190 km s _1 , 440 km s" 1 and —0.75, respectively, in the A-body simulation under the same conditions. 
The agreement is better than it is in ZA (Section 4), probably because ZA does not conserve the momentum (Juszkiewicz 
et al. 1998). This could be a disadvantage of ZA, in addition to the problem of shell crossing. The second-order Eulerian 
perturbation has its own disadvantage that the formulation is too complicated to yield the radial- velocity PDF. Thus ZA and 
the Eulerian perturbation are complementary. Juszkiewicz et al. (1998) showed that {(iiy — (fy)) 3 } approximately scales as 
( v \\)(( v \\ ~ (^y)) 2 )- The negative skewness is induced by the negative average that reflects gravitational clustering of galaxies 
as discussed in Section 5. 

The Lagrangian and Eulerian perturbations are useless in the strongly nonlinear regime, £(r) ^> 1 and accordingly r <C 10 
Mpc at D = 1, for which A-body simulations have shown that the radial-velocity PDF has an exponential tail (Section 4). 
We favor the model of Sheth (1996) and Diaferio & Geller (1996). The exponential tails were reproduced by averaging relative 
velocities of galaxy pairs belonging to the same galaxy clusters. These clusters were assumed to be virialized as well as 
isothermal and have Gaussian velocity distributions with various dispersions according to the theory of Press & Schechter 
(1974). The PDF is symmetric, i.e., (y«) = 0. This is because there was assumed to exist no gravitational clustering. Since 
the A-body simulation yields (i)y) < for r = 1 and 5 Mpc (Section 4), we consider that gravitational clustering still exists 
in practice at these separations. 

7 CONCLUSION 

The radial pairwise velocities of galaxies exhibit a negative average as well as negatively skewed pronounced tails in the quasi- 
nonlinear regime. To understand their origin, we have used ZA, an approximation for gravitational clustering, and analytically 
studied the radial-velocity PDF. 

We have estimated for the first time the separation range where ZA is applicable. For the cold dark matter model of the 
Einstein-de Sitter universe with Ho = 50 km s Mpc , the applicable range is r 3> 10 Mpc at D = 1 (Section 3). Any 
discussion on ZA has to be restricted within this range. 

We have compared our ZA calculation with the result of the TV-body simulation done under the same initial condition 
by Fisher et al. (1994). ZA reproduces successfully the negative average and negative skewness of the radial- velocity PDF 
(Section 4). 

The negative average and negative skewness originate in galaxy clustering according to the initial gravitational potential, 
i.e., the infall vu < dominates over the expansion vu > 0, as pointed out by Juszkiewicz et al. (1998). This information is 
contained in the dependence of the initial velocity dispersion on the initial separation (Section 5). 

The radial-velocity PDF in A-body simulations has a pronounced tail in the negative side. Since ZA cannot fully reproduce 
this property, we have attributed it to gravitational acceleration of clustering galaxies that is not fully incorporated in ZA 
(Section 4). To confirm this conclusion, we would require higher-order approximations, e.g., post-ZA and post-post-ZA. They 
use not only the initial velocity Vi but also the initial acceleration Vi and so on (Bouchet et al. 1995; see also Yoshisato et 
al. 1998). The approximations Pade-post-ZA and Pade-post-post-ZA developed by Matsubara, Yoshisato & Morikawa (1998) 
would be also useful. 
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APPENDIX A: TIME EVOLUTION OF PROBABILITY DISTRIBUTION FUNCTION 

Let us expand the joint PDF, P(v'n , v' ± \r' , t + dt), around the time t using ZA (equations and 0: 



P(ypv' x \r',t + di)Sv\\ 2tyv' ± Sv' ± Anr' 2 Sr' 



P(v»,v±\r,t) + 



BP D 



dP , OP , 
+— vudt+ — dt 
or at 



— «ll + — } dt + 

D ' 



D 



dP ( D 



dv± 



Vj_V\\ 



D 



dt 



(Al) 



1 + 3 — dt I 8v\\2nv±8v±'iivr 2 Sr. 
D ' ~ 



From comparison between equations an d I IAH . we obtain the time-evolution equation of the joint PDF as 



dP 
~dt 



D 
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dP 



dP 



W " dv\\ ^ V± dvi 



+ 3P - 



dP dP 

Vl - a v \\~5 — 

ovu av± 



dP 
dr 



(A2) 



This equation is equivalent to a collisionless Boltzmann equation because ZA ignores encounters among the individual galaxies. 
The first, second and third terms represent the acceleration, rotation and translation of the galaxies, respectively. However, 
(dm ,v±,r) do not constitute the generalized coordinates. This fact yields the third term 3P in the parenthesis of the first term 
(see equation 13 . 

The integration of equation (|A2fl by dv^2irv±dv± yields the evolution of the two-point correlation function £(r,t), 



5(i + {M)) = ¥=4s[Ai + {M))(i)] 



which corresponds to a continuity equation. This is a reasonable result because the galaxy number is conserved as in equation 
0- 



APPENDIX B: PAIRWISE VELOCITIES IN OTHER COSMOLOGICAL MODELS 

Thus far we have assumed the Einstein-de Sitter universe with a Hubble constant Ho = 50 km s _1 Mpc - . The application to 
other cosmological models is straightforward. We only have to know the values of D and D and substitute them to equations 
(0-(|2jl. The functional forms of D are summarized in Suto (1993). 

We compare ZA with A^-body simulations of Magira et al. (2000, see also Suto et al. 1999). They considered three 
models, SCDM (fi = 1.0, A = 0.0, h = 0.50 and V = 0.50), LCDM (Q = 0.3, A = 0.7, h = 0.70 and F = 0.21) and OCDM 
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Figure Bl. Radial-velocity PDFs for SCDM, LCDM and OCDM at r = 22.8 h' 1 Mpc. Solid lines: ZA calculation. Dashed lines: 
W-body simulation of Magira et al. (2000). Dotted lines: Gaussian PDFs with the same dispersions cry (r) as ZA's PDFs. The abscissa is 
in units of km s . 



(fio = 0.3, Ao = 0.0, h = 0.83 and F = 0.25). Here fio is the density parameter, Ao is the normalized cosmological constant, 
and h is defined as i/o/(100 km s _1 Mpc -1 ). Thus SCDM is the Einstein-de Sitter universe, while LCDM is a flat universe 
with a cosmological constant and OCDM is an open universe. The power spectrum of the initial density fluctuation Viiki) 
was from the cold dark matter model of Bardeen et al. (1986, see also Sugiyama 1995): 

7>i(*0 = BD ^%\l^ q) i 1 + 3 - 89 ^ + ( 16 -^) 2 + ( 5 - 46 ?) 3 + ( 6 - 71 ?) 4 ] " 1/2 . (Bl) 

where q = ki /Th. The normalization factor B was determined by comparing the radial- velocity dispersion at infinity cm (r — > 
oo) in the linear theory of equation 11811 with that in the A^-body simulation (591, 606 and 603 km s -1 for SCDM, LCDM 
and OCDM). Fig. Bl shows the radial-velocity PDF at 22.8 ft" 1 Mpc. This separation is within ZA's applicable range. ZA 
and the A/-body simulation are in satisfactory agreement, except for the tails of the PDF that are pronounced in the A/-body 
simulation. In SCDM, the average, standard deviation and skewness are -66 km s" 1 , 589 km s" 1 and -0.22 in the N-body 
simulation, while they are — 88 km s _1 , 551 km s" 1 and -0.18 in ZA. In LCDM, the average, standard deviation and skewness 
are —153 km s _1 , 616 km s _1 and —0.51 in the A^-body simulation, while they are —155 km s _ , 496 km s _1 and —0.38 in ZA. 
In OCDM, the average, standard deviation and skewness are —127 km s , 618 km s _1 and —0.40 in the N-ho&y simulation, 
while they are —136 km s , 512 km s _1 and —0.28 in ZA. 



